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TECHNICAL MEMORANDUM NO. 1211 


THE CHARACTERISTICS METHOD APPLIED TO STATIONARY 
TWO-DIMENSIONAL AND ROTATIONALLY SYMMETRICAL 

GAS FLOWS* 

By F. Pfeiffer and W, Meyer-Konig 


SUMMARY 


By means of characteristics theory, formulas for the numerical 
treatment of stationary compressible supersonic flows for the two- 
dimensional and rotationally symmetrical cases have been obtained 
from their differential equations. 


INTRODUCTION 


The auxiliary means for the theoretical treatment of stationary 
gas flows at supersonic velocity are the characteristics of the 
partial differential equations governing the motion. The Busemann 
graphic methods 1 for the treatment of such potential flows is based 
on the network of the characteristics \ = const., \i = const, in the 
velocity field (u, v— plane ) and the fact, that at corresponding 
points of the velocity and the flow field (x, y-plane ) corresponding 
characteristics are perpendicular to one another. Guderley 2 has 
extended the Busemann method to two-dimensional and rotationally 


*"Die Charakteri st ikenmethode bei stationaren ebenen und rotations— 
symmetrischen Gasstromungen. " Zentrale fur wissenschaftliches 
Berichtswesen der Luftfahrtforschung des Generalluf t zeugmei sters (ZWB) 
Berlin— Adlershof . Forschungsbericht Nr. 1581 , March 20, 1942. 

1 Busemann, A,: "GaBdynamik.” Handbuch d. Experimentalphysik 

Bd. 4, 1. Teil, 1931, pp. 341-460, particularly p. 421 and the 
following pages. 

^Guderley, G, Die charakterist ikenmethode fur ebene und 
achsensymmetrische lib er schall stromungen . Jahrbuch der Deutschen 
Luftfahrtforschung 1940, pp. I 522 - I 535 . 
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symmetrical, turbulent flows. On the basis of physical considerations, 
he defined the changes dX, dp, which the quantities X and p 
experience- in an advance from one point of the flow field to an 
adjacent point along a characteristic. 

For cases more general than Guderley postulated, it is 
possible to get the differential equations for the characteristics 
by pure mathematics directly from the differential equations for 
the gas motion. They have a very simple form, so that they are 
also suitable for the graphic, numerical treatment of special 
problems. 


SYMBOLS 


x, y rectangular coordinates in the flow plane 
v velocity vector 

u, v rectangular coordinates in the velocity plane 
w, -a polar coordinates in the velocity plane 

p density 

p pressure 

e entropy 

T absolute temperature 

cp,c y specific heat at constant pressure, or at constant volume 

a sonic velocity 

a x critical sonio velocity 

i heat content 

i Q heat content at w = 0 

a Mach angle 

V, u', v*, a', p*, p' dimensionless quantities instead of 

w, u, V, a, p, p 
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ABBREVIATIONS 


a - l xi p f , T 

= m (a’ 2 ), 

q _ X 

s - cvp 

S 2 

= s ft = s* - T , 

iC 

1 

x 

i 

c 

- k = 

R = %, - 

c_ 


n — 

sin a cos 

, JS. - , 

a 

y 



I: THE TWO-DIMENSIONAL FLOW 


From the fundamental hydrodynamic equations for compressible, 
stationary flow 3 


(z^i = - p grad p 
div (pv) = 0 


( 1 ) 


for the velocity vector x, the pressure p, and the density p, the 
following equations are obtained for two-dimensional flow with 
velocity coordinates u, v 






In addition for an ideal gas there occurs the equation 


p = Ce s p k 


( 2 ) 

(3) 

(4) 


where S » s the entropy, k = the ratio of tho specific 

heats, C represents a constant, which is derivable from the 
equation of state 


E 

p 


ss 


E T 


(4a) 


3 Webster, A. G-., and Szego, G.: Partielle Differentialgleichung^n 
der mathematischen Physik, 1930, p. 34. 
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because of the first and second lave of thermodynamics (E * °p - %) 
T = absolute temperature). Lastly there is the equation 



+ V< 


gg 

3y 


o 


( 5 ) 


which expresses the hypothesis that the entropy along each individual 
streamline is constant. From ( 4 ) tar the (variable) sonic velocity • a 


a 2 - o fcCeSpk -1 
op 


k£ = kRT 


further 


dfi 

dx 


a£ + a 2 & 
k ox ox 


\ 

l 



( 6 ) .... 


(6a) 


If equations (2) are written thusly: 

A + v& - + vdn + ga + gs = o 

ox ox ox dy p ox k ox 

u Su + v^-u^ + uj2+a2ba + ai3s =0 
dy dy dy dx p dy k dy 


multiplication by dx or dy and addition (with w 2 3 u 2 + v 2 ) 

vdw + ~ dp ' + ^dS = $idy-v dx) ( g& - ) 

p k loy dx/ 


that is, for laminar flow the following prevails everywhere 

2 2 

w dw + ~ — dp + ~~ dS = 0 
P fc 


( 7 ) 
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for turbulent flow, on account of (5) > on each streamline 

p 

w dv + dp =0 (7a) 

P 


On each streamline p = (Cl constant) and a 2 = 

therefore ~ dp = d^ E. j , Accordingly, it follows from (7&) 

by integration for each streamline 


w 2 

2 


+ 



( 8 ) 


or 


w£ 

2 


+ i 


= ^o 


(8a) 


k a 2 

if I = CpT = -■ — ■ - RT e is the heat content of the ideal gas. 

The constant i Q is, in general, different from streamline to stream- 
line. The condition (8a) also holds in the final form if the stream- 
line passes through a compression shock *; first,, only the region that 
does not Include a compression shock will be considered. 


Equations (2), (3), and (5), where and in (2) are 

ox dy 

replaced in accordance with (6a), together with 




du 


du 


dy 


du « ~ dx + ~ dy, dv = ^ dx + 


dx 


dp 


dy 


dp 


dx 


dS 


dv 

dy 


dS 


dy 




dp = dx + ~- dy, dS = ^ dx + ~- dy 


dx 


dy 


dx 


dy 


(9) 


from eight linear equations for 
with the matrix 


du du dv dv dp dp dS ds 

mmmm 9 S ~~ 

ox dy ox by bx by bx by • 


4 Busemann, A.,' (See footnote 1) p. h-33 . 
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du ?u dv V 5S . as 

dx Hy ^x ■ <$ <§y S 'Sy 


i 



The characteristics are here defined as those curves on the integral 

au cii av av ap 

surfaces sought j along which the derivatives x-, - N ~ , 

■s . ^ n. cx dr c.c cy ox 

~ do not result from, these eight equations as definite, 
dy ox ay 

along which, therefore, discontinuities in these derivatives are 
possible. In order that the quantities 

s 5 ax ay ax. ay ox ay 

can he ambiguous, the determinants of eight roue of this 

ax' oy 


matrix must equal zero. Setting these determinants equal to zero 
furnishes the four systems of characteristics (projection of the 
first counted twice) 




u dy -v dx = 0 

dS = 0 

w dw +■ a 2 — = 0 
P 

(v 2 - a 2 ) dx 2 -2uv dx dy + (u 2 
$idy -vdx) (udv - (vdu) + 

K V 


( streamlines ) V 

-J 

- a 2 ) dy 2 = 0 
+ £j~) a 2 (u.dx +vdy) 



( 10 ) 


(ID 


(Mach waves) 
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As a result of the introduction of polar coordinates with u = w cos -6, 

v = w sin ■d in the u, v.-plane and of the Mach angle a through 

sin a = - there is obtained 
w 


uiv -vdu « dd, u du + v dv = w dw 

and from the first equation (11) 

4JT , -uv £ a) 

a 2 - u2 


For w > a, that is, for the case of supersonic flow, which is the 
present problem, the characteristics become real. On account 'of 

V — s T tan a equations (10) and (11) become 
u dx+ vdy 


^ = tan <9 
dz 

dS ■ 0 

dw + a sin a = 0 


= tan (d + a) 


dfl 


sin a cos a 


dp dS n 

+ t + tt = 0 




(10a) 


(11a) 


(These equations are still valid, if the equation of state is in the 
more general form p * f(p, S) provided in the second equation (11a) 

jj: is replaced by ^ 

* 

The equations are made dimensionless, as 'a result of introducing 
w - w°w*, a = w°a f , p = p°p f , where w° is any reference velocity, 
p° any reference density. After setting in' p' = a and S* = p 

the characteristic equations become 
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' flZ * tan fl- 
dx 

dS * 0 


r 


dw* + a r sin a da = 0 J 


(10b) 


= tan (3 T a) 

y Si + da + dS* 

sin a cos a 



(lib) 


with sin a = =r. 

w* 

In the point-by— point ’numerical determination of the flow 
phase with the aid of the characteristics the Massau method 5 is 
made use of. 

Let all the phase quantities at points 1 and 2 of. the flow 
be known and provided with the subscript 1 or 2. A point 3 of 
the flow with its phase quantities is obtained thusly: According 

to the first equation (lib) the position xj, y 3 of 3 is obtained 
from ?3 - yi = tan (3 - a^) (X3 - x^) and 73 - 7g * tan fa 2 + «2) 
(x^ — Xg) graphically or numerically. From the second equation (lib) 
are obtained these equations as differences 


1 

sin oq_ cos ctj. 


($3 -’§1) + (a 3 + S3 Z ) - + Bf 


and 


± (« 3 -fl a ) + (c 3 + S 3 X ) * a 2 + 8/ 

sin ag cos a g 


From this (numerically) ^3 and 03 + S3 X . By means of the first 
equation (10b) 


7 3 - y m = tan <3 (*3 " atm) 


s Encylclopadie der Math. Wissenschaften II, 3, 1, p. 162. 
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either graphically or numerically point m is reacned, from 
point 3 by going in the direction of the streamline related to 3* 
m being in the connecting line 1-2, whose phase quantities are 
computed by linear interpolation from 1 and 2 . Then S 3 X - and, 

therefore 0 3 also, is obtained. The third equation ( 10 b) gives 

w 3 * - w m * = - a* m sin ctm (03 - 0 m ) 
from which w^* follows. 

With this x 3 , y 3 , -8 3, 03, S3*, W3* are to be had. For further 

calculation a r 3 and 0.3 from 0.3* = kCe^ 3 p 3^— d and sin cc.3 = — 3 - 

are used. From the first approximations obtained in this way, 
better values can be obtained by iteration after the manner of 
Massauj often estimated, possible extrapolated values can be used 
for the first approximations. 

In all cases graphic representations can be given for simpli- 
fication of the calculation, for example, because 


3in a cos a /l - (Jjj® 

the ray system A = constant and a = constant from the origin of 
the a*, w* plane, or the system of straight lines S = const, 
in a plane with log a and log p in the coordinate directions. 


NOTE 


In an isentroplc laminar motion, dS 
therefore, everywhere and 


0 and 


Su dv 
dz ciy 


v dw t 



0 


( 7 b) 


The equations for the characteristics systems ( 11 a) become 
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g - tan ( ti + a) 

T dfl - cot a = 0 


( 12 ) 


If P and P* are two corresponding points in the x, y— and the 
u, v-planes and two pairs of corresponding characteristics with th°ir 
tangent directions are drawn through them, then 


tan 7 = — - — = + cot a 
dw/dv 


also 


7 = 90° + a 

3 =s d.+ 7 = 90 0 + ■d + a 


By the first equation (12) 


S = $ . + a 


In addition, there is the fact that is the hasis for the Busemann 
graphic method for isentropic, laminar motion, that in two corresponding 
points of the u*- v— and the x, y-planes the characteristics are 
mutually perpendicular crosswise. 


II. TBS ROTATCONALLY STMMETRICAl, SLOW 

Equations (2) result from equations (1) again hy application of 
cylindrical coordinates 6 and considering the rotational symmetry 
while the term enters in equation (3) on the left side. In 
addition to that, equations (4) and (5) are again valid. 

In setting up the -characteristics, instead of 0 in the 
ninth column, third place from the top appears now and, as a 
result, the characteristic equations become 


6 Frank, Ph.y and 7. Mises, R . : Die Differential- raid ' „ 

Integralgleichungen der Mechanlk und.Fhysik I, second ed., 1930, p. 86. 
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u dy *- v-xlx = 0 
dS = 0 

¥ dw + ^ = 0 


(v 2 _ a 2 ) fa.2 - 2 u v dx dy + (u 2 - a 2 } dy 2 

' i 

fcdy -7 dx ) 2 ^ + (Sidy -vdx) fcdv-vdu) 

+ (J& + a 2 dx + v dy) = 0 



With the polar coordinates w, -8 in the u, v— plane and the 
Mach angle a there is obtained because of 

JL » « JL (* tan a) 1 L°2$ .± dr ( w cos 8 tan ( 8 + a) - w sin* 

a^udx+vdy^a 2 z V 

k 

= g£ QPP* , ■ 

x cos a cos (8 t - a) 



tan 8 


dS = 0 


dw + a sin a da 


0 


^ = tan (a T a) 
dx 

- £L§ + d£ 

sin a cos a p 


+ 


dS , dx cos 8 

k x cos a cos (8 + c) 



(13a) 


providing x £ 0. 


III. SPECIAL CASE-: IDENTICAL i 0 FOR ALL STREAMLINES 

Guderley 7 restricts himself to the practically important special 
case that i D in equation- (8) is the same constant for all stream- 
lines, that is, there is a curve (not a streamline) in the region 
of flow, which is intersected by all streamlines and along which i 0 
has a fixed value. 

It follows from equation (8) then 


w dw + a da = 0 

k — 1 


and from this, because of 


a 2 = kCe s pk— 1 


wdw + a 2 t 2 - + ^ — 

P k_- 1 

7 Guderley, (see footnote 2) p. I 524. 


dS = 


0 
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or 


w dw + a 2 ^+^dS = — 


k(k - 1) 


as 


in the entire region of flow. 


s 


With the resultant equation 


dp . as 1 . 2 da 

T~~ STF'~17 F^t “ 


the second equation (11a) becomes, if, in addition, m (a*2 ) = t 

- — Lui — as - as* + aT * o (lie) 

sin a oos a 

As a result, the numerical work is simpler than in the general 

case. 

Again, all phase quantities at points 1 ana 2 (fig. 1) are 
known. The equations 

y 3 ~ yi = tan ( 8 1 - ai) (x 3 - X]_) 

and 


73 - 72 = ta O (tf 2 + Q&) (x 3 - X2) • 


again give (graphically or numerically) the position x 3 , y 3 of 
point 3. From the equations (lie) expressed in differences' 5 


-Ai (^3 - $ i ) - ($3 - a> x ) = 0 


and 


Ag (*3 -*2) - (» 3 -® 2 ) = 0 


®.Croeoo, L. : Eine neue' Strcntfuhktian fur die Erfoschupg' der'Bewgung 
der Case mi t Rotation, ZAMM 17, 1937, pp. 1-7, particularly p. 2, , 


equation (l”)!?. 


p 
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with A = — and S 1 — t => $ becomes 

sin a cos a 


. = H, li_.+.fag g .i , 

3 A X + Aq 


$3 ■= - A ± (d 3 - %) = ® 2 + Aq (*3 - -Sg) 


By means of the streamline through j, m. i3 obtained, aa in th 9 general 
case, and with S 3 s - = sJ 5 , T 3 » and a* 3 . Because of equation ( 8 ) 
the third equation (10b) is not needed this time. Figure k furnishes 
W 3 and the quantities 0,3 and A 3 needed for further calculation . 

For comparison with Guderley the additional term which is to 
be pu* on the left side of the second equation ( 12 ) for the potential 
flow, is determined in order to obtain the left hand side of the 
second equation ( 11 c) 


+ da - sa g S°5 , » as* 


sin a cos a J( jc . 2 sin a cos a da 

+ — m — t 


for the turbulent flow. This additional term is 


a to . fl -j B a P.ff as* + 4a + 00 t a ^ = - sin. a__cos a d$ 

k - 1 k - 1 a w k - 1 k 


cot a + cot a ~ 
w w 


— cot a 


a c 


ds 


w^ (k - ljkoy 


_ ds a 2 

— ■" COt CL — 

w 4 * KR 


— cot a JL ds 8 
w2 


1 2 

since Cy (k — 1 ) = R end « t 

KR 

In the rotationally symmetrical case there is a further additional 
term; according to equation ( 13 a) it is 


•Guderley, (see footnote 2), I 325 , equation (16a) and ( 16 b). 
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dx s ip a cot' d 
x cos (d + aj 


Introducing 


d: 


I dx 2 + dy 2 = 


dx 

cos (i T a) 


it becomes 


di . 

— sin a cos d 


as in Guderley if it is borne in mind that there, as opposed to 

our designations, y is used instead of x and d there is the 
angle with the axis of rotation, 90° — d , therefore, in the case 
of this report. 


It. COMPRESSION SHOCKS 

The treatment of a compression shock appearing In the two- 
dimensional flow with the variables used in this report is given 
for the example treated by Guderley 11 . 

For any point of the shock line let w be the velocity 
(parallel to the x-axis), p the pressure, p the density to the 
left of the shock line, w (coordinates u, v) velocity, p pressure, 
p density, a sonic velocity to the right of the shock line, 
let p 0 , 0 ^ be the pressure and density of the flowing gas at 
rest, a z rhe (constant) critical sonic velocity and let 




u 


i 


t 




o 


* p 1 - 

■y 



a* 

> 


— — c> 

10 Guderley (gee footnote 2), p. I 526. 


3 1 Gudez ley (gee footnote 2), figure 14. 
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then the flow to the left of the shock line' is given by ’v?* J 
and x p f follows from p* = P ,l£ . The following equations 12 
hold for passing through the shock line. 


.i2 = 


= (17* - 





(shock polar; it is the same curve here for all points of the shock 
line) . 


= p r + p* (W* - u f ) w' grjrj 

(15) 

XlB. + _i_ a’ 2 = i k ..± 3- 
2 * F-T 2 k - 1 

t 

(16) 

p » « a «2 p» 

(17) 

k + 1 


p» = e®p lk 

(18) 


(Equation of state* If the entropy is set equal to zero for p Q * p 0 .) 

Point 1: The direction of the velocity after the shock is 
tangent to the profile at 1; as a result* from the shock polars 
are obtained (fig. 6) wi*, u^'* v^*, and the direction of 
the shock line at 1* perpendicular to AB. 

Then with that* p. 1 from (15)# ai* from (16)* pi* from (17)* 

1 a-j * 

from (lo)* and ai from sin ai = 

W-j_ f 

Point 2: Point 2 is taken close to 1 on the shock line and all 

quantities (except X2 and y2) are taken as at point 1 (fig. 7). 

Point 3: The coordinates x^* y? of the point 3 which is the 

intersection of the characteristics 2-3 leaving from 2 with the 
contour of the profile* are obtained (graphically) by means of 


12 Busemanh (See footnote 1) p. 436. 
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- y2 51 ‘tan (8 2 - °^) (*3 ~ * 2 ) 


Since 1—3 is a streamline, 83* * S-j 2 $ 3 is the direction of 
the tangent to the contour at 3, Tg is computed from 

- Aq (*3 -*2) - (S 3 X - S 2 X ) + (t 3 -t a ) =0 

with t 3> a’g is obtained from 2n (a* 2 ) =r, w f g from 

u j2 . |S 

Sal + fa 1 iuli 

2 k - 1 2 k — 1 

and ag with it. 

Point 4 : With w f and w ! 2 the direction of the shock line 

at 2 is established (it still agrees with the direction- at 1 ) . The 
direction of the characteristics out. from 3 gives, by means of 


n - y 3 a tan (a 3 + a 3 ) - Xg) 


the point 4 (x^, y^) (graphically) as the intersection of this 
characteristic with the direction of the shock line: d 1. is 


estimated and w *4, u’k, v f k 
shock polar; then, with that, 

P'k by ( 17 ), from ( 18 ), 

t 4 from Tjj t ]n (a*^) 2 . 


are obtained graphically from the 
P'k ( 15 ), by (16), 

ak from sin ai. = 

^ ... 4 w*k 


Supplementing thiSj there is the second equation of the 
characteristics from 3: 


Ag («4 - $g) - (S4 X - Sg x ) + (rjj -Tg) =0 

from which a new is computed, with which an iteration is 

carried out. 

Proceeding from point 4 along a characteristic of one family 
to a farther point on the contour, the quantities at this point 
are computed analogous to the manner at 3, from here proceeding 
along, a characteristic of the other family to a farther point of 
the shock line and there computing the quantities analogous to the 
manner at 4 . 
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As soon as the steps become too large, it becomes necessary to 
interpose points] new points in the region between shock line and 
contour are computed by the method in III. The related quantities 
for the points of departure a, b are determined by interpolation 
between 3 and 4 (fig. 8 ). 


Translated by Dave Fe Ingold 
National Advisory Committee 
for Aeronautics 
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Figure 3. 




Figure 4. 






